Simulating Hamiltonian dynamics using many-qudit Hamiltonians and local unitary 

control 



O 

o 

>>■ 



O 
(N 



> 

in 



in 
o 

o 

Oh! 
I 

> 

a • 
^ : 



Michael J. Bremner,!- Dave Bacon,2.3'|t| and Michael A. Nielseni-'i'^^Itl 

'School of Physical Sciences, The University of Queensland, Queensland 4 072, Australia 
^Institute for Quantum Information, California Institute of Technology, Pasadena CA 91125, USA 
'"^ Department of Physics, California Institute of Technology, Pasadena CA 91125, 
^School of Information Technology and Electrical Engineering, 
The University of Queensland, Queensland 4012, Australia 



USA 



When can a quantum system of finite dimension be used to simulate another quantum system 
of finite dimension? What restricts the capacity of one system to simulate another? In this paper 
we complete the program of studying what simulations can be done with entangling many-qudit 
Hamiltonians and local unitary control. By entangling we mean that every qudit is coupled to every 
other qudit, at least indirectly. We demonstrate that the only class of finite-dimensional entangling 
Hamiltonians that aren't universal for simulation is the class of entangling Hamiltonians on qubits 
whose Pauli operator expansion contains only terms coupling an odd number of systems, as identified 
by Bremner et. al. [Phys. Rev. A, 69, 012313 (2004)]. We show that in all other cases entangling 
many-qudit Hamiltonians are universal for simulation. 



PACS numbers: 03.67.-a, 03.65.-w 



I. INTRODUCTION 



A. Overview 

One remarkable aspect of Nature is that it can be mod- 
eled by equations whose solution may be obtained by al- 
gorithmic means. This empirically observed fact allows 
us to construct physical theories that make predictions as 
to how Nature will behave. Of course, while we can sim- 
ulate Nature, our capacity to do so is limited by the way 
we choose to perform the simulation and the complex- 
ity of the system to be simulated. Feynman's landmark 
paper on quantum computation Jlj] discussed the appar- 
ent inability of classical computers to efficiently simulate 
quantum systems and suggested that a quantum com- 
puter might succeed where classical computers fail. In 
this paper we study a class of simulation protocols moti- 
vated by the example of quantum computation. In par- 
ticular, we examine the following question: given a com- 
posite system with a finite-dimensional Hamiltonian and 
the ability to perform arbitrary local unitary operations, 
what other Hamiltonians can we simulate? 

The simulation of quantum systems by quantum com- 
puters is a topic that has attracted considerable atten- 
tion. A considerable literature (see and references 
therein) addresses the question of how to adequately sim- 
ulate physically interesting closed quantum systems. Is- 
sues of particular interest include the complexity of pro- 
tocols for simulating initial states, simulating evolutions, 
and for extracting physically important information from 
the final state of the computer. Each of these issues must 



be addressed in any comparative study of quantum and 
classical computers, and their capacity to simulate Na- 
ture. 

While state preparation and measurement are vital 
elements of any simulation of a quantum system, we 
focus in this paper on the simulation of evolutions of 
systems. Hamiltonian simulation protocols using single- 
qudit unitary operations as an additional resource have 
received considerable attention in recent years due to 
their relationship with various models of quantum com- 
putation. One of the more noteworthy advances was the 
discovery that all two-body Hamiltonians can simulate 
all other Hamiltonians on the set of qudits that they en- 
tangle, when combined with single-qudit unitary opera- 
tions UlllllSlllSlillillltll- This body of work 
also demonstrated that these Hamiltonians could be used 
to efficiently simulate any other two-body Hamiltonian 
that acts on the network of qudits they entangle. This in- 
cludes a Hamiltonian that can implement the CNOT oper- 
ation, thus implying that all entangling two-body Hamil- 
tonians and single-qubit unitary operations are universal 
for quantum computation. 

The results and tools used to study two-body Hamilto- 
nian simulation have been applied fruitfully to several re- 
lated problems. There is now a considerable literature on 
time-optimal strategies for simulating two-qubit Hamil- 
tonians and quantum gates; see, for example. [^ll3 . IT^ 
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references therein. This body of investigation has led to 
interest in applying these theoretical results to practical 
proposals for quantum computation '2^. 

More recently, studies have focused on using systems 
with many-qubit interactions for Hamiltonian simulation 
and gate-synthesis |23, 11^ HO, |^ ■ A number of these 
papers have investigated the structure of systems with 
many-body interactions for the pur poses of gate synthe- 
sis and algorithm design |2j |23, yfl ISJ ■ Several authors 
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have recently examined the effects of many-body inter- 
actions in quantum dot [s^l and optical lattice |33. l3^ 
systems. 

For the purposes of this paper we are most concerned 
with the work in (2Sl|, where the authors established 
which Hamiltonians with many-qubit interactions are 
universal when combined with the ability to perform 
arbitrary single-qubit unitary operations. In a similar 
vein we examine which Hamiltonians with many- qudit 
interactions are universal when combined with arbitrary 
single- qudit unitary operations. Our final result is a strik- 
ing generalization of the conclusion in . showed 
that the only class of non-universal entangling Hamil- 
tonians on qubits are the odd entangling Hamiltonians, 
i.e., those Hamiltonians whose Pauli operator expansion 
contains only terms coupling an odd number of qubits. 
Furthermore, showed that the odd entangling Hamil- 
tonians can all simulate one another, so there is a sense in 
which there are only two types of many-qubit entangling 
Hamiltonian. Remarkably, in this paper we will see that 
when the systems involved are not all qubits, this struc- 
ture actually simplifies, with all entangling Hamiltonians 
capable of simulating all other entangling Hamiltonians, 
i.e., we show that apart from the many-qubit case, there 
is only one type of many-body entangling Hamiltonian. 

Our primary concern in this paper is with questions 
of universality in many-qudit systems, without regard to 
the issue of complexity. Thus, when we say a set of re- 
sources is universal on a set of qudits, we are stating that 
these resources can be used to simulate any Hamiltonian 
on those qudits, without implication that this simulation 
is efficient or inefficient. This is in contrast to the notion 
of universality for quantum computation which requires 
that any universal set of resources can simulate a stan- 
dard gate set with a polynomial overhead in the number 
of qubits used. That said, it is often possible to exploit 
the structure of certain classes of many-body Hamilto- 
nians to develop efficient simulation algorithms. For in- 
stance much headway can be made in developing efficient 
Hamiltonian simulation protocols that use /c-local Hamil- 
tonians by adapting the methods developed for systems 
of qubits in j23| to systems of qudits. 



B. Terminology and statement of results 

Before turning to the discussion and proof of the main 
results of this paper it is helpful to introduce some ter- 
minology. Generally, we will use the term qudit to de- 
scribe any quantum system with a finite-dimensional 
state space. As an example of our usage, a three- 
qudit system might contain a two-dimensional system (a 
qubit), a five-dimensional system, and a four-dimensional 
system. 

We are interested in the properties of the Hamilto- 
nian dynamics of an n-qudit system. As we will see, a 
great deal can be said about the properties of a Hamil- 
tonian simply by examining its structure in a suitable 



representation. In |29j, the authors found that the uni- 
versality properties of a many-body Hamiltonian acting 
on qubits could be identified by expanding it in the Pauli- 
operator basis, i.e., tensor products of X, Y, Z and /. 
In this paper we expand upon this analysis by examin- 
ing the properties of an 71-qudit Hamiltonian written in 
a d-dimensional generalization of the Pauli basis. 

An arbitrary Hamiltonian on n qudits can be uniquely 
written as 

H = J2ho.H^, (1) 

a 

where ha are real coefficients and each Ha in the expan- 
sion is a tensor product of Hermitian operators acting on 
the individual quits, 

n 

Ha^lS)Hi^\ (2) 

where Ha^ acts on qudit j, and is either the identity 
operator, or one of a set of traceless Hermitian matri- 
ces known as the Gell-Mann matrices. The Gell-Mann 
matrices generalize the Pauli matrices, and thus this ex- 
pansion is a generalization of the expansion for qubits 
used in '29*1 . The Gell-Mann matrices for a d-dimensional 
quantum system consist of: (a) d — 1 matrices of the form 

W„ = , / (j2 \b){b\ ~ (m- l)|m)(m| I , 

y'm{m- 1) J 

(3) 

where 2 < m < d; and (b) the Pauli- like matrices: 

Xat = -^i\a){b\ + \b){a\) (4) 

Y,, = _l(|a)(6|-|6)(a|) (5) 

where 1 < a < 6 < d. These act as the Pauli X and Y 
on the two-dimensional subspace spanned by the vectors 
I a) and We sometimes refer to the Wm matrices as 
Cartan subalgebra elements of the Gell-Mann matrices, 
since they span a Cartan subalgebra of the Lie algebra 
su{d) generated by the Gell-Mann matrices. However, 
it is worth emphasizing that we do not use any special 
properties of Cartan subalgebras, and the reader does not 
need to be familiar with the properties of Cartan subal- 
gebras to follow the details of the paper; our use of the 
term is a convenience of nomenclature only. Note that 
the Gell-Mann matrices are traceless and Hermitian, and 
form a complete basis for traceless Hermitian matrices. 

The representation Eq. ^ is useful as it highlights 
which qudits interact and which do not. In particular, 
given a term Ha let Sa be the set of qudits upon which 
Ha acts non-trivially, that is, the set of qudits for which 
Ha is traceless. We say that the qudits in Sa are coupled 
by Ha and refer to Ha as a coupling term. We also say 
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that Ha is entangling on the set Sa- More generally, we 
say that a Hamiltonian H is entangling on some set of qu- 
dits if it is not possible to partition this set of qudits into 
two non-trivial sets S and S such that every term Ha in 
the expansion of H couples either a subset of S or of S. In 
graph-theoretic language, if the qudits corresponded to 
vertices on a hypergraph and the couplings corresponded 
to hyperedges, the condition that the Hamiltonian is en- 
tangling on a set of qudits is simply that the hypergraph 
connects this set. As such we say that a Hamiltonian 
connects the set of qudits it entangles. 

Our strategy for demonstrating universality in this pa- 
per is to show that some set of resources is capable of 
simulating another set already known to be universal. In 
particular, we make reference to two theorems that cat- 
egorize large classes of Hamiltonians as universal up to 
single-qudit unitary operations. The first was mentioned 
in the introduction: two-body entangling Hamiltonians 
are universal for quantum computation 10]. Using the 
terminology just introduced, this theorem may be stated 
as follows y\Sj\: 

Theorem 1. Suppose H is a two-body Hamiltonian, that 
is, every coupling term in the Gell-Mann expansion of H 
couples at most two qudits. If H is entangling on a set of 
n qudits (that is, the coupling terms in H connect these 
qudits) then evolutions of H together with single-qudit 
unitary operations are universal for quantum computa- 
tion on these n qudits. 

The second universality theorem that we use involves 
Hamiltonians acting on sets of qubits that have coupling 
terms that may couple more than two qubits. This the- 
orem is stated |2^ : 

Theorem 2. Suppose H is an arhitary entangling 
Hamiltonian on a set of n qubits. Evolutions of H and 
single-qubit unitary operations are universal on those n 
qubits if and only if the Gell-Mann (i.e., Pauli) expan- 
sion of H contains at least one coupling term that couples 
an even number of qubits. 

Theorem |2 tells us that for a Hamiltonian acting on 
qubits alone to be universal, it must have a coupling term 
acting on an even number of qubits. If H does not con- 
tain such a coupling term then we shall call it an odd 
Hamiltonian, since all its terms couple an odd number of 
qubits. What can the odd Hamiltonians simulate? This 
question was also answered in [Igj : 

Theorem 3. Let H be an odd entangling Hamiltonian, 
that is, every term in the Gell-Mann expansion of H cou- 
ples an odd number of qubits. Then H and single-qubit 
unitaries can simulate any other odd Hamiltonian on the 
n qubits. 

pol l also demonstrated that the Lie algebra generated 
by the odd entangling Hamiltonians on n qubits (and lo- 
cal unitaries) corresponds to the Lie algebras so(2") and 
sp(2"), for even or odd n respectively. Furthermore, p9l | 



showed that the odd Hamiltonians can be made universal 
with appropriate encodings. 

In this paper we demonstrate that if a Hamiltonian is 
entangling on a set of qudits, then this Hamiltonian is 
universal on those qubits, when assisted by local unitary 
operations. The only exception to this result is the spe- 
cial case when the Hamiltonian is an odd Hamiltonian 
acting on qubits only. 

C. Outline 

Theorem n shows that if a Hamiltonian connects a set 
of qudits with two-qudit couplings then this Hamilto- 
nian is universal with single-qudit unitaries. Our strategy 
in this paper is to show that a many-body Hamiltonian 
(that isn't one of the odd qubit-only Hamiltonians) con- 
necting a set of n qudits can simulate a two-body Hamil- 
tonian connecting the same set of qudits. This is done 
by defining a series of simulation protocols, each identi- 
fying broad classes of Hamiltonians that any entangling 
Hamiltonian can simulate, until we arrive at the eventual 
result. 

The structure of the paper is as follows. In Section UTI 
we introduce some simple general simulation techniques 
that are used often in this paper. Section HiH introduces 
a simulation technique known as term isolation. This 
simulation technique allows us to simulate any particu- 
lar coupling term, Ha, that is present in the Gell-Mann 
expansion of H, thus isolating the term. In Section IIVI 
we show that given some term coupling k qudits, we can 
simulate new coupling terms that couple fewer than k qu- 
dits. We also discuss the limitations on this type of sim- 
ulation. Section Ivl examines how we can use a term that 
couples k qudits to simulate a coupling between two qu- 
dits. Finally we prove the main result of the paper: that 
the only non-universal class of entangling Hamiltonians 
is the class of odd Hamiltonians. This is argued through 
an exhaustive demonstration that all n-qudit entangling 
Hamiltonians other than the odd many-qubit Hamiltoni- 
ans are indeed universal. 



II. SIMPLE SIMULATIONS 

In this section, we review some simple Hamiltonian 
simulation techniques studied in previous papers H, ^ |M 
H 0, H i, ES m [13, and that will form the basis for 
our later results. By a Hamiltonian simulation we mean 
a sequence of evolutions due to our system Hamiltonian, 
H , which is assumed fixed, interleaved with single-qudit 
unitary operations. The goal is to approximate (to arbi- 
trary accuracy) evolution according to some other Hamil- 
tonian. If that is possible for some desired Hamiltonian 
we say that Hamiltonian can be simulated. The theory of 
Lie algebras and Lie groups ensures that the techniques 
decribed in this section exhaust the set of possible sim- 
ulations that can be performed given some Hamiltonian 
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and single-qudit unitaries. 



A. Conjugation by a unitary operator 

A quantum system with Hamiltonian H evolves in time 
via the unitary operation e~*^*. Say we are also given 
the ability to perform some unitary operation, U ^ and 
its inverse, C/^ . Then performing the sequence of unitary 
operations Ue~'^^*W = g-^UHuU^ that we can 

simulate an evolution according to the conjugated Hamil- 
tonian UHW. In this paper, as we have given ourselves 
the ability to perform arbitrary single-qudit unitaries, we 
will often conjugate a Hamiltonian by unitaries of the 
form U = Ul(g)U2<E)■■■<E)Un■ 



li. Simulating linear combinations 

Suppose we can simulate two different Hamiltonians, 
Hi and H2. Then we can simulate the sum of these 
Hamihonians, since g-^^i'^e"'-^"^ « g-i(ffi+i/2)A 
small A, and with successive evolutions we can simu- 
late the Hamiltonian Hi + H2 for an arbitrary time t. 
Imagine that we could evolve our system by a whole set 
of Hamiltonians, H, and their negatives^. It follows that 
we can simulate arbitrary linear combinations of any of 
the elements of H. 



C. Simulating commutators of Hamiltonians 

Another simple simulation protocol that can be 
performed is the simulation of a commutator of 
two different Hamiltonians. This is possible as 

simulate Hi , H2 and their negations we can simulate the 
commutator of these Hamiltonians. 



D. Simulating Hamiltonians that couple the same 
qudits 

Consider the general expression for a Hamiltonian act- 
ing on a system of qudits in Equations (^) and |(2Jl, and 
recall that Ha couples a set of qudits Sa ■ We now intro- 
duce a theorem from |3| to show it is possible to use Ha 
and single-qudit unitaries to exactly simulate any other 
coupling term that couples the set of qudits 5*0,: 



Given that we can simulate H, it turns out always to be possi- 
ble to simulate —H, using single-qudit unitary operations. This 
follows from Equation 1111 . later in the paper, which shows how 
to express —H as a sum of terms of the form UHU'^ , where U 
are local unitary operations. By the methods of simulation we've 
already introduced, it follows that —H can be simulated. 



Theorem 4. Let A and B he any two traceless Hermi- 
tian operators in d dimensions and assume that 0. 
There is an algorithm to find a set of at most unitary 
operators, [/„, and constants Cn > such that: 

A = ^c„f/„Si7t. (6) 

n 

Key to proving this theorem is a result from the theory 
of operator majorization, Uhlmann's theorem 35] . Al- 
though we do not need the theory of majorization in this 
paper, for the benefit of readers familiar with majoriza- 
tion, we make the following summary remarks. Recall 
that Uhlmann's theorem tells us that if P ^ Q (that is, 
P is majorized by Q) then P = Y^^PnUnQUl^^ for some 
unitary operators [/„ and some Pn that form a probabil- 
ity distribution. The proof of Theorem^jin follows by 
showing that A -< cB for some positive constant c. 

Any coupling term Ha in is a tensor product of 
traceless terms acting on Sa- If we replace B in Theo- 
rem 0] by the individual tensor factors appearing in Ha, 
then we see that we can simulate any A that is a tensor 
product of traceless Hermitian operators acting on the 
same set Sa- This result will be extremely useful in the 
remainder of this paper. It tells us that if we can simulate 
some coupling Ha , we can simulate every other coupling 
on the same set of qudits. 

III. TERM ISOLATION 

In Section III Dl we saw that any coupling term. Ha, 
in the expansion H = J^a^aHa (Equation (^), could 
be used to simulate any other coupling term that entan- 
gles the same set of qudits. If we have a Hamiltonian 
that is simply a coupling term on a given set of qudits, 
we can immediately say a great deal about what can be 
simulated with that Hamiltonian. In general we do not 
have this luxury of interpretation. Instead, some general 
Hamiltonian, H ~ "^a^a, has many different coupling 
terms that couple many different sets of qudits. Term 
isolation is a simulation technique that uses H and single- 
qudit unitaries to simulate any particular term Ha in the 
expansion of H alone. 

Term isolation allows us to think about in a different 
way, showing that the ability to simulate H is equivalent 
to the ability to simulate the coupling terms {Ha} in- 
dividually. Thus, we can perform our analysis entirely 
in terms of the set {Ha} and still encapsulate all of the 
Hamiltonian simulation properties of H . Given that the 
elements of the set {Ha} have a much simpler structure 
than a general H, term isolation is a powerful tool for 
analysis. 

We now show that term isolation can always be per- 
formed. If we demonstrate that we can use H and single- 
qudit unitaries to simulate some Ha coupling an arbitrar- 
ily chosen set of qudits, then we know from Section III Dl 
that it can be used to simulate any other term coupling 
the same qudits. 
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Without loss of generality we may assume that the 
term being isolated is of the form 

k 

where k is the number of qudits in the set Sa ■ To see that 
there is no loss of generality in assuming this form, note 
that we can always relabel the qudits in Sa so that they 
are the first k qudits in the system, and any operators 
Xab or Yab in Ha are equivalent under local unitaries to 
W2. 

Any term in the expansion of Hfj, that isn't the 
term Ha that we wish to keep, is different from Ha in at 
least one of three ways. Either: 

Case 1: H/3 has terms acting non-trivially on qudits out- 
side of Sa, the set of qudits upon which Ha acts. 

Case 2: Hp acts on a strict subset of Sa- 

Case 3: Hp acts on the same qudits as Ha but is a ten- 
sor product of different elements of the GcU-Mann 
basis. That is. Ha ^ Hp, even though Hp couples 
the set Sa- 

Each of these cases identifies a special difference between 
Ha and Hp. In the following sections these differences are 
exploited to define simulations that remove undesirable 
terms. 

As we have previously stated, every simulation in this 
section may be represented as a sequence of linear com- 
binations, commutators and conjugations by local uni- 
taries. We often denote a sequence of operations of this 
type on a Hamiltonian, H, by a scripted letter. For ex- 
ample, in Section llll Al we define the depolarizing chan- 
nel, which is a linear combination of conjugations by lo- 
cal unitaries, and write ~ Hd to symbolize the 
depolarizing channel acting on H, resulting in the simu- 
lated Hamiltonian Ho . The action of I? on defines a 
simulation. We can also compose simulation techniques, 
so, for example, in Section IIII Bl we define a simulation 
T[Hd] = Ht- 

A. Case 1 

We begin by noting the identity 

J2 UpJUl = d tT{J)I, (8) 

where J is an operator acting on some qudit of dimension 
d and the sum is over all elements of the d-dimensional 
Pauli group^, where we omit repeated summation when 



^ The proper ties of the d-dimensional PauH group were extensively 
studied in l36l . We will not use any further sp ecia l properties 
of this group and refer the interested reader to |3fij for further 
information. 



two elements in the Pauli group differ merely by a phase 
factor. We note there is a simple extension of Equa- 
tion for multiple- qudit systems, 

= Dtr{J)I^", (9) 

where the superscripts indicate the different qudit sys- 
tems, of respective dimension d^-'-', D = d'^^\..d'^'^^ is the 
dimension of the combined system, / represents the ap- 
propriate identity operator for each subsystem, and the 
sum is over conjugations by all elements of the Pauli 
group for each qudit, again omitting repeated sums over 
elements that are the same up to a phase factor. 

We define the simulation = Ho to be the 

multiple-qudit depolarizing channel acting on the n — k 
qudits that aren't coupled by Ha , 

= Hd- (10) 

Ha acts on the first k qudits of an 71-qudit system, that 
is, the set Sa- If we examine the simulated Hamilto- 
nian, Ho, we find from Equation ^ that any terms Hp 
in H that act non-trivially on qudits outside the set Sa 
are eliminated. The simulation leaves the coupling term 
Ha unchanged except for an unimportant positive scal- 
ing factor. Thus we have removed all the Case 1 terms 
Hp from the Hamiltonian, and need only consider the 
remaining Case 2 and Case 3 terms. 

B. Case 2 

The Hamiltonian Ho is a linear combination of terms 
that couple the set of qudits Sa or some subset of Sa- 
lt turns out that we can use another extension of Equa- 
tion ||SJ) to simulate a Hamiltonian, Ht, that only has 
terms that couple the set Sa- In Equation JSJ, if J is a 
traceless operator we find that the right hand side of the 
equation is zero. Noting that / is an element of the Pauli 
group, we find 

^ U,,JUl = - J, (11) 

which always holds for traceless J . Using single-qudit 
unitaries from the Pauli group we consider the following 
summation, 

^ ([/w®c/f )(j«®j(2))(f^a)^t/f )t. (12) 

If J(i) and J(2) are traceless, this expression is equal to 
® If j(^) is traceless and J^^^ is the identity. 
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this expression is equal to — [(d^-^))^ — 1]/ (g) J^'^\ With 
this in mind we define a simulation: 

(13) 

Performing T*^^^ for j — 2, ...,fc, only terms that couple 
the same qudits as Ha are not eliminated. So performing 
the following sequence of simulations, 

t[Hd] = r(^)[r('=-i)[...[r(2)[i/c]]...]] = Ht (w) 

the simulated Hamiltonian, H^, is a linear combination 
of terms that couple the same qudits as Ha- 



C. Case 3 

We have shown how to simulate a Hamiltonian Ht 
that only contains terms which couple the same qudits 
as Ha- To eliminate the remaining terms we define the 
following operators that are both unitary and Hermitian, 

d 

Za^I-2\a){a\=Y,\j){j\-2\a){a\. (15) 

Notice that the Za operators commute with each of 
the Cartan subalgebra elements, Wm, in Equation 
Hence, each of the Za will also commute with Ha as it is 
a tensor product of elements of the Cartan subalgebra. 
Further notice that Za anti-commutes with Xim and Yim 
if a — I or a = m and commutes otherwise. We can 
use this fact to define a simulation that eliminates terms 
with Xirn and Yim operators present in Ht- We define a 
simulation 

Z(J^[H]=H + Z^pHZi^\ (16) 

where the superscript j indicates a Za operator acting on 
the jth qudit, with identities acting elsewhere. If there 
exists any term with an Xim or Yim operator on the jth 
qudit, and such that a = ^ or a = m, then this term 
will be eliminated from Ht by the simulation Za^HT]- 
Expanding on this idea we can eliminate every term on 
the jth qudit that has the form Xim or Yim by performing 
the following simulation: 

^ Z^,'\zII^^[...[z['\Ht]]...]] (17) 

where d is the dimension of the jth qudit. The effect of 
this simulation on Ha is simply to rescale it. Now, if we 
perform the simulation Z^^' for each qudit in Sa, 

Z[Ht] = zW[Z('=-i)[...[Z(i)[i/T]]...]] = Hz, (18) 



all that remains in the newly simulated Hamiltonian, Hz, 
is a linear combination of terms that commute with the 
Cartan subalgebra elements. We have now simulated a 
Hamiltonian with no X- and K-type terms. 

Hz is a linear combination of terms that are ten- 
sor products of operators from the Cartan subalgebra. 
Consider the unitary representation, P^^\-k), of the per- 
mutation group Sbj-i that permutes the elements of 
the diagonal basis of the Cartan subalgebra, |a), for 
a — 1, . . . , 6j — 1 on the jth qudit. When a > bj we 
have P<-J^TT)WaP'-J^HTT) = Wa- When a < bj, we find 
that the effect of conjugating Wa by a permutation op- 
eration is to shift around the diagonal elements of Wa- 
Now, we can eliminate any terms in Hz that contain an 
operator Wi"'' with a < bj by performing the simulation 

V'-^^[Hz]= J2 P^^^HzP^^^^- (19) 
(i) 

This works because Wa is a diagonal, traceless oper- 
ator and the permutation, V^^\ distributes each of the 

diagonal elements of W^^ equally. The effect of P on 

(i) 

terms Wa acting on the jth qudit and with a > bj is 
to simply scale them by a factor of {bj — 1)!. Performing 
the following simulation, 

V[Hz] = V^''\v'^''-^\-[V^^^[Hz]]-]] = Hp (20) 

we produce a Hamiltonian Hp that is a linear combina- 
tion of terms that couple the same qudits as Ha and are 
tensor products of operators Wa with a > bj. 

In Section III CI we pointed out that it is possible to 
simulate a Hamiltonian proportional to the commutator 
of two Hamiltonians that are both simulatable. Now, we 
note that the commutator —i[Wa''\Xh^^ih^] = if a > 

bj- If a = bj we find -i[wl'\ Xb^^it,] - ^=n,-i6,. 

We can make use of this distinction to find a way to 
remove the unwanted terms from Hp- We define the 
simulation 

(21) 

Then if we perform the following sequence of simulations, 

X[Hp] = A'W[A'('=-i)[...[A'(i)[iJp]]...]] = Hx (22) 

we find that Hx ^ {(^'j=iYb^-i.b,) «> J®""'', up to some 
unimportant but non-zero constant multiple. We have 
now simulated a single coupling term that couples the 
same qudits as Ha- Recall in Section Hi Dl we noted that 
a coupling term can be used with single-qudit unitaries to 
simulate any other term coupling the same set of qudits. 
So, we can use Hx and single-qudit unitaries to simulate 
Ha, the desired term. Thus we have demonstrated that 
it is possible to isolate Ha from H- 
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IV. SIMULATING NEW COUPLING TERMS 

Term isolation shows that the abihty to simulate a 
Hamiltonian H = haHa is equivalent to the ability 
to simulate the set of coupling Hamiltonians, {-ffa}, given 
single-qudit unitary operations. Additionally, we learnt 
in Section III Dl that given Ha and singlc-qudit unitaries 
we can simulate any coupling term that couples the same 
qudits as Ha- So far we have not presented any way of 
simulating some coupling term that couples a different 
set of qudits than any of the terms in the set {Ha}- In 
this section we will take a key step towards a proof of uni- 
versality, showing how to use single-qudit unitaries and 
a term Ha coupling k qudits in order to simulate a term 
that couples fc — 1 qudits. 



A. Evaluation of commutators 

In [2^ it was shown that if Ha coupled qubits, its 
capacity to simulate other coupling terms depended on 
the number of qubits that it coupled. More specifically, 
it was shown that if Ha coupled k qubits and k was 
an odd number, then Ha couldn't be used with single- 
qubit unitaries to simulate a coupling term that coupled 
k — 1 qubits. One way of seeing why this is true is to 
examine the commutator of two Hamiltonians, [-ffa, Hfj]^ 
that couple the same set of qubits Sa - It is easy to show 
that the commutator \Ha , Hp] 7^ if and only if there 
are an odd number of locations in Sa where Ha and Hfj 
differ. From this restriction it is possible to prove, as was 
done in , that coupling terms coupling an odd number 
of qubits can only ever simulate other Hamiltonians that 
have odd couplings. 

What is different when not all the systems are qubits? 
The purpose of this subsection is to investigate the com- 
mutator of two specially chosen couplings Ha and Hp 
that couple the same set of qudits, Sa- In the case of 
qubits, it is not difficult to convince oneself that when 
Sa contains an even number of qubits, the commutator 
[Ha, Hp] is either zero, or else couples a set of qubits 
that is a strict subset of the original set Sa- We will 
show by an explicit calculation that when one or more 
of the systems is not a qubit, it is possible to choose Ha 
and Hp so that the commutator [Ha , Hp] contains terms 
coupling the entire set Sa - Remarkably, we will see in the 
remainder of the paper that this is the key fact that sim- 
plifies the study of universality when not all the systems 
are qubits. 

We begin by choosing Ha = ^ab ^^'^ ~ 

^'^^iX^i^l where for all j we set h ^ h' - (We assume 
initially that all systems are of dimension 3 or greater.) 
Given these forms for Ha and Hp , what does [Ha , Hp] 



look like? We find 

-<^^^{xiil^^Y,^^^ (23) 

This expression contains Hermitian and skew-Hermitian 
terms. Upon expansion of the above expression we find 
that all of the Hermitian terms sum to zero, leaving only 
a sum of skew-Hermitian terms remaining. These terms 
correspond to a sum of tensor product terms containing 
odd numbers of lf,b' terms. All of the terms couple the 
entire set Sa ■ It is easy to verify that this sum is always 
non-zero, simply by inspection of the coefficients of the 
relevant terms. 

So far we have only considered the case where we could 
choose to simulate Ha and Hp for xj^J and , b ^ b' - 
We can only do this when each subsystem has dimension 
d > 2. If we have subsystems where d — 2, the situation 
changes slightly, but the results are similar, provided not 
all of the subsystems are qubits. 

For every j where the qudit has dimension d > 2 we 
choose ni'^ = X^^f^ and i?^^' = with h ^ b'. For 
every j where the qudit has dimension d — 2, we choose 
hL^^ = X, and iJ^^' = Y. Provided Ha and Hp do 
not couple qubits exclusively, a straightforward calcula- 
tion along lines similar to that already done shows that 
[Ha, Hp] is a non-zero sum of terms, each of which is 
skew-Hermitian and couples all k qudits. The only sub- 
tlety in the calculation is the need to analyse separately 
the cases where there are an even number of qubits in 
the set Sa, which gives rise to a commutator which is a 
non-zero sum of tensor product terms containing an odd 
number of Ybb' terms, and the case where there are an 
odd number of qubits in the set Sa, which gives rise to 
a commutator which is a non-zero sum of tensor product 
terms containing an even number of Yi,b' terms. 

B. Simulating identity operators 

Given some term, Ha , coupling a set of qudits Sa , we 
show how the results on commutators just obtained allow 
us to simulate other coupling term that couples a subset 
of Sa with just one qudit removed. More precisely: 

Lemma 1. Given the ability to evolve via Ha = 
{^"^-^ Ha \ which couples k qudits, and local unitary op- 
erations, it is possible to simulate H' such that 

H' = I®H^, (24) 

provided Ha does not couple qubits exclusively. The cou- 
pling term H^ may couple any k — 1 qudit subset of Sa , 
subject to the constraint that the subset not be qubits ex- 
clusively- 
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Proof: Given Ha we can simulate any other coupling 
term, Hp = (^^^iH^ , that acts non-trivially on the 
same set of k qudits, Sa- We label the qudits so that 
Sa consists of qudits 1, . . . , fc, and so that our goal is to 
simulate a coupling on qudits 2,. . . ,k, i.e., the goal is 
to remove qudit 1. To this end, we choose H^^'' so that 



^(1) 



Hj^\ Note that, by assumption, the set 2, 



does not contain qubits exclusively. Evaluating the com- 
mutator, we find: 



i[Ha,Hf3]=i{Hi'yf. 



\H, 



j=2 



(25) 



Setting N = 0;^2 Hi'\N' = 0;^2 and applying 

Equation JSJl to the first qudit, we see that it is possible 
to simulate 



H' = il® [TV, N'] 



(26) 



Finally, we note that as N and N' don't act exclusively 
on qubits, our earlier results on commutators show that 
we can ensure that [N, N'] is a non-zero linear combina- 
tion of terms that couple Sa^ less the first qudit. Term 
isolation allows us to simulate one of the coupling terms 
in [N,N'] alone, i.e., H" — I ® H^, as required. □ 



V. UNIVERSALITY 

Theorem ^ stated that if a set of qudits is connected 
by a Hamiltonian, H, with two-body interactions, then 
evolutions by H and single-qudit unitaries form a univer- 
sal set of operations on that set of qudits ^ ^^t of 
2-qudit coupling terms connecting the same set of qudits 
is also universal as they can simulate a two-body Hamil- 
tonian on the set of qudits. We prove in this section the 
main result of this paper: that a generic Hamiltonian, H , 
entangling a set of qudits can simulate a set of 2-qudit 
coupling terms connecting the qudits, and is thus univer- 
sal. The only exception to this rule is the case where H 
is a sum of odd coupling terms, as discussed in and 
summarized in Theorems |21 and 13 in the present paper. 

We begin by proving Theorem |S1 which shows that a 
coupling term. Ha, that couples a set of k qudits, Sa, can 
be used to simulate a set of 2-qudit couplings that con- 
nect the set Sa- This implies that Ha and single-qudit 
unitaries are a universal set on the qudits Sa- We con- 
clude with Theorem El showing that an arbitrary entan- 
gling Hamiltonian on n qudits is universal for the qudits 
it entangles. 



A. Theorem 151 Using a term coupling many qudits 
to simulate a term coupling two qudits. 

Theorem 5. Suppose Ha — ® j=i ^a'"' couples k qu- 
dits. Then Ha and single-qudit unitary operations can be 



used to simulate a set of two-qudit couplings connecting 
every qudit coupled by Ha, provided Ha does not couple 
qubits exclusively, and k > 1. Thus Ha and single-qudit 
unitaries are universal on the set of qudits coupled by Ha - 

Proof: Without loss of generality we may label the 
systems so that Ha couples systems 1 through k, and 
system 1 is not a qubit. Fix j in the range 2 through 
k. Applying Lemma ^ repeatedly, we see that we can 
simulate a Hamiltonian coupling system 1 and system j . 
It follows that Ha and single-qudit unitaries are universal 
on the set of qudits coupled by Ha- □ 



B. Theorem 151 Which Hamiltonians are universal? 

With Theorem[Slin mind, we now prove that the only 
non-universal set of entangling Hamiltonians is the set of 
odd Hamiltonians acting on qubits alone. 

Theorem 6. Single-qudit unitary operations, and evo- 
lutions via a Hamiltonian, H , that connects a set of n 
qudits, is a universal set of operations on those n qudits 
if and only if H is not an odd Hamiltonian acting on 
qubits alone. 

Proof: The forward implication follows from Theo- 
rem |21 as does the reverse implication when all systems 
are qubits. Thus, all that needs proof is the reverse impli- 
cation in the case when H is an entangling Hamiltonian 
that does not act exclusively on qubits. We will show how 
to construct a set of two-body couplings that connect all 
n qudits. 

To construct this set, begin by picking a system that 
is not a qubit, and label it system 1. We will explain 
how to construct a set, S, of systems to which 1 can be 
coupled via a two-body interaction. We begin by setting 
S — {1}, and aim to add in other systems that can be 
coupled to 1 via two-body interactions. Our strategy is 
to show that provided S is not yet maximal, i.e., does 
not yet contain all n qudits, then it is always possible to 
add an extra qudit into S. 

To see this, suppose S is not yet maximal. Then it is 
always possible to pick a qudit j inside S and a qudit k 
outside of S such that H contains a coupling term Hjk 
which couples systems j and k. (Other systems may also 
be coupled by Hjk-) In the case when either j or k is 
not a qubit. Theorem [S] shows that a term coupling just 
j and k may be simulated. Theorem ^ implies that we 
can also simulate a term coupling system 1 and fc, and 
so system fc may be added to S. 

The other possible case is when j and k are both qubits. 
In this case, suppose without loss of generality that Hjk 
has the form X^^^^ (g) X^*^^ (g) . . ., where the superscripts 
label the systems. We may also simulate the coupling 
since system j is in S. Taking the commuta- 
tor of these two couplings, we see that we may simulate 
couplings of the form ® Y^^'' ® X^'^^ (g) . . .. Applying 
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Theorem 13 we see that it is possible to simulate a two- 
body coupling between system 1 and k, and thus system 
k may be added to S. □ 

VI. CONCLUSION 

We have demonstrated that many-qudit Hamiltonians 
combined with local unitary operations are always uni- 
versal for simulation on any connected set of subsystems 
upon which the interactions act nontrivially, provided 
that Hamiltonian is not an odd Hamiltonian acting on 
qubits. This result is rather intriguing and elegant, espe- 
cially in the light of the general lack of broad results for 
many-body (as opposed to two-body) problems in quan- 
tum information science. In the study of pure state bi- 
partite entangled states, for example, a single unit of 
currency, the maximally entangled state, has been iden- 
tified and the fungible nature of this currency has been 
established. On the other hand, a similar currency and 
set of fungible transformations has not been identified 
for systems consisting of more than two parties. Given 
this difficulty in understanding the structure of quantum 
states, it is quite remarkable that, with the exception of 
odd entangling Hamiltonians, all of the different many- 
qudit interactions are equivalent. Even in the case of 
odd entangling Hamiltonians, universal simulation can 
be achieved using an encoding which wastes only a single 



extra qubit of space [29j. Thus there is a real sense in 
which, for simulation, all interactions have been created 
equal. 

Part of the simplicity of our result stems from our focus 
on universality for simulation as opposed to universality 
for quantum computation, which requires that issues of 
efficiency be taken into account. When one adds the re- 
quirement of efficiency of simulation, then problems of 
universality become much more difficult: indeed this is 
perhaps one of the fundamental problems in the study 
of the computational complexity of quantum circuits. A 
well-developed theory of efficient simulation is a task of 
great importance and, judging from the difficulties en- 
countered in proving lower bounds for problems in classi- 
cal circuit complexity, this task is probably an immensely 
difficult problem. This paper can be seen, however, as a 
necessary precursor to any attempt to advance this pro- 
gram. 
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